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Abstract 

The new excluded volume hadron gas model by Singh et al. is critically discussed. We demon- 
strate that in this model the results obtained from relations between thermodynamical quantities 
disagree with the corresponding results obtained by statistical ensemble averaging. Thus, the model 
does not satisfy the requirements of thermodynamical consistency. 
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I. INTRODUCTION 



In order to take into account the effect of particle repulsion at short distances the excluded 
volume procedure was proposed by van der Waals in 1873. In a system with a fixed number of 
particles N this procedure corresponds to a substitution of the system volume V by the available 
volume V — vN, where v is a proper particle volume. This proposal was further supported by 
the calculations within statistical mechanics for a system of hard spheres of radius r at small 
particle number densities. In this case, v equals the particle hard-core volume 47rr 3 /3 multiplied 
by a factor of four. The excluded volume procedure leads to the equation of state: 



p(V - vN) = NT , (i; 



or equivalently 



p (1 — vn) = nT , (2) 



instead of the ideal gas relation 



Pid = n id T . (3) 



In the above equations, p and T are pressure and temperature, respectively, and n = N/V is 
particle number density. 

The hadron e cLS IS cL system with conserved charges (baryonic number, electric charge, 
strangeness), but with a variable number of particles. Thus, for the hadron gas the excluded 
volume procedure ([1]) with a variable number of particles is needed. The first step is to de- 
velop the excluded volume procedure in the grand canonical ensemble (GCE), as in the GCE 
a summation over number of particles is performed. 1 This step appeared to be not trivial. 
The correct excluded volume procedure in the system with a variable number of particles was 



first introduced in Ref. (8[. Several other proposals |9l-lllj do not satisfy the thermodynamical 
relations. 



Recently, the new excluded volume model has been developed by Singh et al. 



13] (we refer 



to it as NEVM). In the grand canonical partition function the authors substitute the system 



1 The excluded volume procedure with a variable number of particles is also needed for the hadron gas in the 
canonical and micro-canonical ensembles. For these ensembles the conserved charges are exactly fixed, but 
the number of particles still remains a variable quantity. 
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volume V by the available volume V(l — Ylj v j n j)-> where sum is taken over different particle 
species, rij is the j-th particle number density, and Vj is the proper volume parameter of the 
j-th hadron. Contrary to the statements in Refs. [l-7] this procedure is not thermodynamically 
consistent. In order to prove this we consider the simplest example of the system with a single 
particle type obeying Boltzmann statistics. 



II. THERMODYNAMIC AL CONSISTENCY 

As a first step let us discuss general requirements of thermodynamical consistency. In GCE 
the partition function is given as the sum over states i: 



Z(V,T,„) = Y>p(^^) 



(4) 



where iVj and are the number of particles and system energy in z-state, T and \i are the 
system temperature and chemical potential, respectively. The system pressure is calculated as 

p = L i n Z(y,T,») , (5) 

whereas particle number density n and energy density e are found using the thermodynamical 
relations 

dp 

n{T,fi) = , (6) 

Thus, thermodynamical average quantities for the number of particles and system energy are 
equal to: 

N = Vn(T,fi) , (8) 
E = Ve(T,fi) . (9) 

In order to fulfil the thermodynamical relations (jHJ) and ((Zj) a correct structure of the partition 
function as a sum over the system states i is required. The quantities T and fi may only enter 
Z(V,T,fi) in the very definite way, i.e., only as given by Eq. (TJJ. Under these conditions the 



3 



GCE statistical averages 



(iV> 
(E) 



Z(V, 



z(y,T,fjL) 

can be calculated using T- and /^-derivatives as 

z{v,T,») ^ ^ ex H r 



cxp 



(JV) 
(E) 



z(y,T it i) 







J2 ( t2 ^ + T i'— 



T 



(10) 

(11) 

(12) 



T 



T^ — + Tfi— Z(V, T, fj,) = V 



dp dp 



(13) 



dT ^dfx 

in agreement with thermodynamical relations ([6]) and (I7j), i.e., the thermodynamical averages 
N (jSJ) and (jUJ), and statistical averages (N) (ITU1) and (i?) f lTTj) . are equal to each other: 



(N) = N 



(E) 



DO 



(14) 



III. NEW EXCLUDED VOLUME MODEL pHZfl 

NEVM [l-7 



7j assumes the following expression for the GCE partition function 

fiN\ (V - vN) 



Z<y,T,n) = E-p(^) Z(V,T, N) = £exp(^) 

N=0 ^ ' N=0 ^ ' 



„iV 



AT! 



exp [exp(/i/T) (V - vN) z] 



where 



z(T) 
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k 2 dk exp 



(jfe2 + m 2)l/2 



T 



(15) 



(16) 



is the so-called one-particle partition function. In Eq. ([16]) . m and g are the particle mass 
and the degeneracy factor, respectively. The quantity A" in Eq. (115p is an average number of 
particles calculated with Eq. (jBJ). From Eq. ( TT5l) one finds the pressure: 
T 



p = — In Z(V,T,fi) = Texp(/j/T) (1 - oti) z = (1 - vn) p id 



(17) 



Note that 



n id (T,fi) = exp(p/T)z(T) 



(18) 
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is the particle number density in the ideal gas (i.e., at v — 0), and the ideal gas pressure pid is 
given by Eq. ([3]). The particle number density n(T, /i) is calculated in NEVM [1KZ| using the 
thermodynamical relation ([6]), which together with Eq. f lT7|) yields: 

n = I 1 — vn — vT— ) n i( i . (19) 



dfx J 

The pressure ( IT7I) and particle number density ( [19]) in NEVM are connected by the ther- 
modynamical relation ([6]). Note, however, that in NEVM relation ( 1T71) which connects p and 
n is different from that given by Eq. ([2]). The thermodynamical relations fl6]) and (JTj) do not 
correspond in NEVM to correct statistical averaging of the number of particles ( ITU]) and system 
energy (lUJ. The canonical partition function in ( [15]) reads: 

with A" given by Eq. (jSJ). A presence of A/", which is a function of T and //, in the partition 
function ffT5]) destroys the correct structure (j3]) of the GCE partition function, thus, it leads to 
a violation of the thermodynamical consistency. The thermodynamical relations ([6j and (JTJ) do 
not correspond to correct statistical averaging in NEVM, i.e., 

(AO (E) ±E. (21) 

Therefore, the thermodynamical relations and ([7]) used in NEVM to calculate particle 
number density and energy density are in contradiction with statistical averages (fTU]) and (|lip . 
This is because a presence of A^ = Vn(T,fi) in the partition function ([15]) . which leads to the 
redundant terms in Eqs. ([12]) and ( [T3]) . These incorrect terms come from the derivatives of A^ 
in respect to \i in ( [12]) and in respect to T in ( [13]) . 



IV. EXCLUDED VOLUME HADRON GAS MODEL [8] 
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is based on the equation for the grand canonical partition 



9(V - vN) z N . (22) 



The excluded volume model 
function: 

'liN\ (V - vN) N 

3xp y 

N=0 

The form of the available volume (V — vN) with the theta function 6(V — vN) instead of V — vN 
in (|15p causes some technical problems for the summation in ([22]) over N. On the other hand, 



AT—n \ / 



just this form of the available volume makes the formulation thermodynamically consistent. It 
is easy to see that the thermodynamical relations (EJ) and (171) are in agreement with statistical 
averaging ffTUl) and ffTTj) for the partition function ff22l) . In calculations one uses the Laplace 
transform of (J22j): 

U V CJqj^ V S ) Zj \ V , ± , /JL) — 7 

'0 



/*oo 00 

T, /i) = / rfV exp(y S ) Z(V, T, /i) = V 



AT=0 



exp(— 1> s) z exp(/x/T) 



at=o 

~1 AT 



^(^-wAT) 7V exp(-s1/) 



s - exp(- us) n id (T,fjL) 



(23) 



The system pressure in the thermodynamic limit can be then found from the farthest-right 
singularity s* of Z (J23"j) as a function of s: 



T 



(24) 



The singularity s* is the pole singularity for (|2"3"]) . which leads to the transcendental equation 
for the pressure: 



p = exp (-vp/T) p id . 

Using the thermodynamical relation one finds 

n id exp (-vp/T) 



n 



(25) 



(26) 



1 + vn id exp(-vp/T) 

Note that Eqs. (EMZED are different from Eqs. (ED and of NEVM pJ-tZj| . Combining 
Eqs. fl25H26p and using the relation p id = Tn id one obtains after simple algebra which coincides 
with Eq. 02]). 



V. CLOSING REMARKS 



There are great difficulties in a formulation of statistical theory of hadrons from the first 
principles. Models in this field are of a phenomenological origin. An important constrain in 
model formulation is the requirement of thermodynamical consistency. The thermodynamical 
consistency can be violated in different ways. Two examples are illustrated by considering the 
excluded volume model. In the first example, the inconsistency appears when one substitutes 
the system volume by the available volume in all ideal gas thermodynamical functions. This 
was done, e.g., in Refs. js-ll]. The available volume is a function of temperature and chemical 
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potential, thus, the redundant terms appear, and the thermodynamical relations ( ISTTl) are not 
fulfilled after this substitution. The second example of the inconsistency can be seen in the 
new excluded volume model of Refs. lRTj. In this model one substitutes the system volume 
by the available volume in the GCE partition function. Then the available volume is used to 
calculate the pressure function only. Other thermodynamical functions are calculated using the 
thermodynamical relations (161171) . However, the available volume being a function of T and /x, 
destroys the correct structure of the GCE partition function. This leads to a violation of the 
thermodynamical consistency, i.e., results obtained from thermodynamical relations (ISTTl) and 
from statistical averaging are different. 

Several comments concerning the excluded volume model [8| should be added. The (effective) 
proper volume v used in this model is a model parameter. For a "rigid" balls it is a ball volume 
multiplied by a factor of 4. For "soft" hadrons this factor will have a value between 1 and 4. 
However, in most cases the proper volume v is treated as a free parameter fitted to data. 

A simplest case of a gas of identical particles was considered in the above presentation for 
didactic reasons. The consistent formulation of the excluded volume hadron gas model can 
be easily extended to many particle species with different proper volume parameters (see, for 
example, Ref. 

The effect of relativistic contraction is expected to be small in the hadron gas. This is 
because the hadron gas temperature is significantly smaller than almost all hadron masses and 
thus hadrons are non-relativistic. This is not true for pions, their masses are comparable to 
temperature. But a predominant majority of pions comes from resonance decays, and resonance 
masses are again significantly larger than the gas temperature. Most important relativistic effect 
in the hadron gas model is a variable number of hadrons. But this feature is exactly the main 
point of the consistent formulation ( 12211261) of the excluded volume hadron gas model. 
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